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The characterization of the "most connected" nodes in static or slowly evolving complex networks 
has helped in understanding and predicting the behavior of social, biological, and technological net- 
worked systems, including their robustness against failures, vulnerability to deliberate attacks, and 
diffusion properties. However, recent empirical research of large dynamic networks (characterized 
by irregular connections that evolve rapidly) has demonstrated that there is little continuity in de- 
gree centrality of nodes over time, even when their degree distributions follow a power law. This 
unexpected dynamic centrality suggests that the connections in these systems are not driven by pref- 
erential attachment or other known mechanisms. We present a novel approach to explain real-world 
dynamic networks and qualitatively reproduce these dynamic centrality phenomena. This approach 
is based on a dynamic preferential attachment mechanism, which exhibits a sharp transition from 
a base pure random walk scheme. 

PACS numbers: 89.75.Hc, 89. 75. Da, 89.65.-s, 05.40.Fb 



I. INTRODUCTION 

The modern study of social networks largely presup- 
poses static or slowly evolving complex networks, as 
clearly indicated by an extensive review of the pertinent 
literature [IHS] ■ In the social and behavioral sciences, this 
static perspective of networks has generated a wide spec- 
trum of formal representations and tools for measuring 
structural and locational properties of social networks 
{e.g., centrality and prestige), social network roles and 
positions {e.g., structural equivalence), and local (sub- 
graph) structures. Statistical models of static networks 
are prevalent Q] H] • The study of networks has been de- 
veloped beyond the social sciences by a substantial num- 
ber of papers that have extended the study of universal 
properties in physical systems to complex networks in 
social, biological, and technological systems; these recent 
advances have generated an extensive research effort in 
understanding the effect of static connectivity patterns 
on various dynamical processes occurring on top of net- 
works 7J. One of the most important concepts in com- 
plex network studies is network centrality, which charac- 
terizes the "most connected" nodes in static complex net- 
works [I] |2]- Among other things, the concept of network 
centrality has led to the study of "scale- free networks" [5] , 
where the degree centrality of nodes is distributed accord- 
ing to a power-law or other long-tail distribution, imply- 
ing the existence of highly-connected nodes called hubs 
[3] [5J. The leading explanation for this uneven distri- 
bution of connectivity, present in many real networks, 
is based on growth and preferential attachment mecha- 
nisms, in which the structural changes of real networks 
are governed by the dynamical evolution of the system 
(see [8] and its many variants [9])- The long-tail charac- 
terization of complex networks has helped in understand- 
ing and predicting the behavior of social, biological, and 



technological networked systems, including their robust- 
ness against failures E] , vulnerability to deliberate 
attacks [101 H2], and diffusion properties [T3]. This 
static network centrality perspective has also led the dis- 
cussion of strategies for immunizing particular individu- 
als to reduce disease propagation [TIHTT] . advertising to 
"opinion leaders" [18], and targeting central individuals 
in terrorist networks [19] , 

The analysis of evolving social networks has long been 
regarded as a "Holy Grail" of network scientists [U [2] [20] . 
A number of models, mostly theoretical with limited em- 
pirical backing, have been considered for analyzing lon- 
gitudinal social network dataj2l]]22 , with networks that 
change either at discrete time points [2"3H2l)] or contin- 
uously over time [2"?ll29| . The complex network com- 
munity has also been interested in the analysis of these 
dynamic networks; their research has mainly examined 
various mechanisms of network formation {e.g., "prefer- 
ential attachment" and its variants [31(9]), and the evolu- 
tion of aggregate properties of the network at the macro- 
scopic and mesoscopic levels {e.g., network size, average 
degree, small-world characteristics, or community struc- 
ture [301 [31]). 

While some of the above studies allow their networks 
to change, in most cases the evolution is assumed to hap- 
pen more slowly than any dynamical processes occurring 
on the network. Real network connections, however, are 
often more fluid than this; even when there exists an 
underlying fixed topological structure, connections be- 
tween nodes can adaptively become active or inactive 
over time [201 . Consideration of the dynamical inter- 
play between the state and the topology of the network 
is pertinent to a wide variety of network types. Exam- 
ples include travel on a transportation network, mes- 
sage exchanges over communication networks, interac- 
tions between molecules that can bind to each other, and 
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the switching of genes "on" and "off": all are networks 
with an underlying structure that can change quickly 
and adaptively with respect to the state of the net- 
work, and whose changes are relevant to the behavior 
and functionality of the system. The generality and 
broad implications of adaptive and dynamical wirings 
in real complex networks has been leading to a rapidly 
growing study and much interest from social scientists, 
physicists, and biologists (see reviews of recent advances 
in [32, 33J). This mostly theoretical research has revealed 
a number of new approaches and insights including self- 
organization of adaptive networks contact pro- 
cesses and epidemiology on adaptive networks [3"?H3"5] . 
and social games on adaptive networks [20] [40], [41] . 

Despite these important theoretical advances, there 
have been very few empirical studies related to the dy- 
namics of large-scale network connections, or theoreti- 
cal models that are validated and informed by empiri- 
cal results [12T|44| . Clearly, the study of data collected 
from real-world time-dependent networks, accompanied 
by theoretical models which address and reproduce the 
properties of said networks, will provide new, and per- 
haps fundamental insights into the mechanisms of dy- 
namic networks. For example, the experiments detailed 
in Ref. [42] detected a dramatic time dependence in net- 
work centrality and the role of nodes, something that is 
not apparent from static (i.e., time-aggregated) analysis 
of node connectivity and network topology. These exper- 
iments studied the time-dependent structure of a large- 
scale email network of over fifty-seven thousand users, 
based on data gathered over one hundred thirteen days 
from log files maintained by the email server at a large 
university. They found, among other things, that the 
daily networks (as well as the total aggregate network) 
were scale-free, but the "hubs" (that is, the most well- 
connected nodes) of these networks changed from day to 
day. In other words, a node which is among the most 
connected on one day, might only have a couple links 
the following day (popular today, anonymous tomorrow) , 
and may not even be among the hubs of the aggregate 
network. This phenomenon, called dynamic centrality, is 
not accounted for by the standard explanation of hub for- 
mation: normally, hubs in networks develop in a system 
where popular nodes tend to attract more links (a process 
called preferential attachment), so that when one partic- 
ular node becomes slightly more popular than the rest, it 
attracts slightly more links from other nodes, making it 
more popular, and so forth in a process of positive feed- 
back. This process typically works over a long period of 
time, as news of a node's new "popularity" is allowed to 
disseminate through the system before the next connec- 
tion is made. In Ref. [42], however, hubs develop over 
the course of a single day, and vanish just as quickly: too 
fast for such a process to occur. Thus, the implications 
of dynamic centrality are in sharp contrast to previous 
complex network research, which depend on such prefer- 
ential attachment-type mechanisms to create power-law 
distributions in empirical complex networks. This re- 



sult draws attention to the inherent dynamic nature of 
adaptive networks which has long been neglected in the 
literature [35] • Similar results [35] were found in the in- 
teractions caused by the spatial proximity of personal 
Bluetooth wireless devices, recording the interactions be- 
tween pairs of students over the period of 31 days. Re- 
sults that are consistent with dynamic centrality behavior 
have also been seen in the social dynamics of free and 
open-source (FLOSS) development projects [15], peer 
prestige in academic hiring networks [46) . trading net- 
works [47] . protein-interaction networks [48] , and tran- 
scriptional networks [IS]- The availability of complete 
(i.e., non-sampled) data that tracks the temporal dy- 
namics of the interactions within the network will allow 
further study of dynamic centrality. 

The existence of dynamic centrality and a dynamic in- 
teraction structure has strong implications for existing 
work related to the attack and error tolerance of com- 
plex networks as well as in network transport. Research 
in epidemics has suggested that an effective disease or 
computer-virus prevention strategy would be to iden- 
tify and vaccinate the high-degree nodes of a network, 
which would inhibit the spread of infection [T4l[T7] . Sim- 
ilarly, "popular influencer" marketing techniques (closely 
related to word-of-mouth or viral marketing) are based 
on the premise that focusing marketing activities on the 
hubs of social networks increases the likelihood of a cas- 
cading adoption of products or services, forming a type 
of social epidemic |18j . The underlying assumption of 
both models is that the network topology is basically 
static: a node which is well-connected or popular now 
will continue to be a hub later in time. In a dynamic net- 
work, however, well-connected nodes can quickly become 
only weakly connected (or even disconnected) tomorrow, 
which invalidates such strategies. Thus, the existence of 
dynamic centrality calls for a radical rethinking of the 
interplay between link/structural dynamics and the dy- 
namical processes underlying the time-dependent com- 
plex networks, and may call for modified strategies in 
predicting and preventing (or encouraging) epidemics. 

Our paper contributes to the understanding of adap- 
tive and dynamic networks by demonstrating a new ap- 
proach to explain and reproduce real world dynamic cen- 
trality phenomena. More specifically, we will investigate 
the type of network seen in Ref. [351 S3] , which we will 
call dynamic scale-free (DSF) networks; these are a series 
of scale-free networks with a scale-free aggregate and dy- 
namic centrality. We will show that DSF networks arise 
from a series of vertex-reinforced random walks on a fi- 
nite scale-free network (called the underlay), governed by 
a "limited" preferential-attachment reinforcement rule — 
limited in the sense that, when making a connection, a 
node only takes into account the popularity of its neigh- 
bors in the underlay. 

The study of random walks on graphs, and reinforced 
random walks in particular, has a long fascinating his- 
tory in its own right J5UJ [ST]. In the context of social 
and complex networks, the nature of random walks and 
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diffusion over small-worid and scale-free networks is also 
a much investigated topic [52Tf59] . Random walks have 
been used in analyzing search and navigation on net- 
works [50] , calculating the "betweenness centrality" of 
nodes in a network |61j , extracting community structure 
from a network [61 j . and sampling a large network as a 
means of estimating its topological properties [521 155] . 
The close interplay between the degree of a node and the 
proportion of time visited by a uniform random walk is 
a well-known result [50] , and has also been examined for 
various real world networks [51] . This body of literature 
largely focuses on the behavior of non-reinforced random 
walks on static network topologies. In contrast to dy- 
namics on networks, a different approach for modeling 
the dynamics of networks, which uses reinforced random 
walks or other stochastic processes that incorporate re- 
inforcement mechanisms, has been considered [20] [22]. 
Most pertinent to our study are models of dynamic net- 
works where individuals choose whom to interact with 
depending on how "popular" their choice is [2U1 [2H] • De- 
spite these theoretical advances, application to real world 
dynamic networks has been limited to parameter estima- 
tion with very small networks that are observed at several 
time points [22] , Here we broaden the scope of this re- 
search in several distinctive and novel ways. First, to the 
best of our knowledge, the use of a reinforced random 
walk has not yet been used as a mechanism of dynamical 
network connectivity that is informed and supported by 
actual high-frequency measurements of large-scale social 
network link dynamics. This shift in perspective leads to 
new insights including evidence that non-reinforced ran- 
dom walk schemes cannot serve as a plausible explana- 
tion for experimental data of dynamic networks, in con- 
trast to reinforced random walk schemes. Second, ran- 
dom walk mechanisms (or other stochastic models) have 
not been used so far to study real-world dynamic central- 
ity phenomena. Finally, new analysis methods, including 
various "ranking" and "centrality overlap" measures, are 
incorporated to facilitate and enable effective evaluation 
of models against large-scale data. 



II. MODEL 

In our approach we start with the aggregate network, 
which is called the underlay; it is reasonable to assume 
that the aggregate network represents the long-term pat- 
tern of contact between individuals, which is often well 
approximated by a power-law distribution [3J. Under- 
lays are generated using the Barabasi-Albert algorithm 
[5] with parameters too — 5 (initial number of nodes; also 
the minimum degree of every node in the underlay) and 
N = 1,000,000 (the size of the network). We then gen- 
erate a series of DSF networks using a vertex-reinforced 
random walk upon this underlay network: starting from 
a randomly chosen node, we add weight (according to a 
scheme described below) to the vertices every time they 
are visited, constantly altering the jump probabilities so 



that the walk will tend to revisit vertices already visited 
[20] [51] . This proposal is inspired by the way messages 
on actual communication networks are part of larger- 
scale conversations. For instance, suppose Alice sends 
an email to Bob, who then responds to Alice and also 
forwards her message to Claire. Claire could in turn ei- 
ther forward it to someone new, reply to either Alice or 
Bob, or do nothing. Assuming these actions are roughly 
equal in likelihood, there is a greater probability that 
Claire sends an email to either Alice or Bob, than to 
some other specific individual; thus, people who are al- 
ready part of the email chain are more likely to become 
part of the chain again. Here, each "daily" network is a 
subnetwork of the underlay, built on top of the underlay 
via the node-reinforced random walk mechanism. 

The specific algorithm is as follows: starting with a 
randomly chosen node, we consider that node's neighbors 
i in the underlay, each neighbor being weighted according 
to 



Wi = l + CV h (1) 



where Vi is the number of visits the node has received 
so far, and C is a parameter of the model. Neighbor i is 
then chosen as a target with probability H; — Wi/ 53 • Wj. 
Once a target is chosen, the link between the current 
node and the target is added to the sub-network, the al- 
gorithm steps to the target, and the process continues 
until the chain reaches the prescribed number of steps 
S. The parameter C captures the tendency of the ran- 
dom walk to revisit vertices — that is, the strength of 
dynamic preferential attachment. When C = 0, one has 
a pure (i.e. unreinforced) random walk, a case we will 
refer to as a control, to demonstrate the importance of 
preferential attachment to our results. It should be noted 
that the vertex-reinforced random walk is known [51j to 
have significantly different behavior from the pure ran- 
dom walk. For example, it is well-known that in a pure 
random walk the probability of visiting a particular node 
becomes proportional to the degree of the node, after 
log N steps. This is not so for vertex-reinforced random 
walks which can become "trapped" within a particular 
subnetwork [ST]: nodes within the subnetwork are vis- 
ited frequently while nodes outside the subnetwork, even 
high-degree nodes, remain unvisited. This effect is criti- 
cal to this model's ability to recreate dynamic centrality, 
but from a practical point of view it can be too powerful, 
trapping the walk in a very small subnetwork (even one 
of only two points). To mollify this effect somewhat, we 
modify the vertex-reinforced rule to prohibit links from 
being traversed two times in a row: nodes and links can 
be revisited multiple times during the simulation, but a 
walk may not jump immediately backwards to the node 
it came from. 
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FIG. 1. A log-log plot of the cumulative degree distributions 
of daily networks, averaged over 25 days, built with single 
chains of 1 million steps: each histogram shows the proba- 
bility that a member node of the daily network has degree 
k or larger. All curves show power-law behavior; the inset 
shows the difference between the power-law exponent of each 
daily network's histogram, and the exponent of its underlay 
(approximately -1.9, though a different underlay was used for 
each value of C). 



III. RESULTS & DISCUSSION 

We will now show that the above reinforcement scheme 
is able to create scale-free networks with dynamic central- 
ity and a scale- free aggregate {i.e., the underlay, scale- 
free by construction), reproducing the main qualitative 
features of DSF networks seen in [32 H3]- We begin 
by constructing daily subnetworks with S = 1 million 
steps atop underlays of N = 1 million nodes. Figure [T] 
shows the cumulative probability degree distribution (av- 
eraged over 25 daily runs) for the cases of no (C = 0), 
weak (C = 1), or strong (C = 10) reinforcement. All 
three distributions show power-law behavior. However, 
the networks with reinforcement (C > 0) have smaller 
slopes and longer tails, which suggests a more heteroge- 
neous degree distribution and a more pronounced scale- 
free behavior. Furthermore, if we fit the histograms to a 
power-law Ax~ a , the exponent a for C = is very close 
to the power-law exponent (a u « 2) of the underlay [65], 
suggesting that the random-walk network is merely tak- 
ing on the scale-free nature of the underlay. The inset 
shows this more clearly: considering 25 daily networks 
constructed for each value of C, we fit the cumulative his- 
togram to a power-law, found the exponent ay for each, 
and subtracted the exponent a u of that network's un- 
derlay. (In this case, we used three different underlays, 
one for each value of C.) We see that the C = curve 
decays at the same rate as the underlay (or even a little 
more quickly); the difference is ad~a u = 0.06±0.02. By 
contrast, the C > cases have ay — a u — —0.76 ± 0.01 
(for C = 1) and -0.73 ± 0.04 (for C = 10). With re- 
inforcement, the daily exponents are smaller by about 
0.75, resulting in a slower rate of decay and relatively 



more nodes of high degree. Note also that these results 
are relatively consistent from day to day. 

Having shown that the subnetworks created by our 
model are scale-free subnetworks of a scale-free aggre- 
gate, we now demonstrate that they exhibit dynamic cen- 
trality. Following 42, 43], we define the "Top-n overlap" 
for two subnetworks of a given underlay as follows: af- 
ter ranking the nodes according to their degree in each 
subnetwork, we count the number of nodes which appear 
among the top n nodes in both networks, dividing the 
count by n to express it as a fraction (or percentage). 
If the overlap is close to 100%, then the two subnet- 
works' hubs are largely identical; thus we expect that 
daily subnetworks exhibiting dynamic centrality (where 
hubs change from day to day) have a relatively low Top- 
n overlap, both with the underlay and with each other. 
(The exact method and formula for calculating the over- 
lap, which deals with the issue of ties for nth place, is 
described in the Appendix.) 
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FIG. 2. The Top-2% overlap, as defined in the text, as a 
function of the number of steps S in the daily network, for 
various values of C. The overlap for C > appears to increase 
at a much slower rate, suggesting that the highly connected 
nodes in the underlay are only modestly represented in the 
daily networks. 



We first used this "top n overlap" to compare indi- 
vidual daily networks with their underlays. In Fig. [2] 
we see that the overlap of the top 2% of nodes (i.e., 
the 20,000 nodes with the highest degree) is dramatically 
smaller with the introduction of dynamic preferential at- 
tachment, once the walk extends further than roughly 
200,000 steps, and the overlap decreases monotonically 
as the preference weight C increases. While the overlap 
of the C = case grows steadily as the random walk fully 
explores the underlay, the overlap for C > increases at 
a much slower rate, suggesting that the highly connected 
nodes in the underlay are only modestly represented in 
the daily networks. 

In Fig. [3j we look at the distribution of the Top-2% 
overlap between the underlay and 10 different daily sub- 
networks, generated with walks of S = 1 million steps, 
and for C — 0, 1, and 10; we also look at the Top- 
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FIG. 3. The distributions of the Top-2% overlap between 10 
daily networks and their common underlay (the peaks labelled 
"underlay"), and between the 45 pairs of those 10 networks 
(labelled "cross"). The bar graphs were constructed by bin- 
ning, with a bin size of 0.01%. The results suggest that 'local' 
hubs vary greatly from day to day. 



2% overlap for the 45 pairs one can make among the 
10 daily networks for each value of C. All the distribu- 
tions are sharply peaked, showing that the overlap does 
not vary much from day to day, and the overlap is again 
greatly decreased by the introduction of reinforcement, 
even from day to day. This demonstrates that the set 
of hubs of the daily network varies greatly from day to 
day, a result consistent with empirical dynamic scale-free 
networks 1321 1351. 



50% 



40% - 



30% 



20% 




10% 

2% 5% 



30% 



FIG. 4. The importance of the number of unique visited 
nodes. The inset shows the fraction of nodes in the underlay 
that are included in a daily network generated with different 
numbers of steps S and with C = 0,1,2,5, and 10. Higher 
percentages of unique visited nodes correspond monotonically 
to lower values of C. The main graph shows the same curves 
plotted versus the Top-2% overlap instead of the number of 
steps. The kink in the C = curve is a finite-size effect. 



An explanation for dynamic centrality in the C > 
case can be found by examining the number of unique 
nodes which are visited by the random walk, which is 



equivalent to the number of nodes in the resulting daily 
network, as is done in Fig. [4j For the pure random walk 
(C = 0), the underlay network (due to its small- world 
nature) is rapidly covered as the number of steps S in- 
creases, as is seen in the inset. The Top-2% overlap 
(shown in the main portion of the figure) is closely corre- 
lated to the number of nodes in the daily network, so as 
the daily network grows, it increasingly adopts the under- 
lay's hubs as its own, ruining dynamic centrality because 
the hubs remain constant from day to day. With re- 
inforcement (C > 0), however, the random walk is more 
likely to revisit nodes it has visited before, and so it takes 
much longer to explore the underlay. This slow rate of 
exploration appears to be the primary reason that the 
overlap is suppressed in Fig. [2j in fact, Fig. [4] shows that 
if we compare daily networks with the same number of 
nodes, rather than the same number of steps, reinforce- 
ment actually results in a slightly higher overlap. 
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FIG. 5. Comparison of the underlay network with daily net- 
works. We show the average rank of a node in 10 daily subnet- 
works (with 1 million steps) versus their rank in the underlay 
network; nodes with higher degree have lower rank, and the 
top-ranked node has rank 0. In the case of ties, all tied nodes 
are assigned the rank that is the average value of the range 
they cover. The inset shows, not the average, but the best 
daily rank that node has over the 10 days. The C > case 
is remarkably egalitarian - every node has an opportunity to 
become a hub in some daily network. 



We have so far looked at the behavior of the top 20,000 
nodes in the underlay and their modest role in the daily 
networks, but we can also compare the rankings of all 
nodes in the underlay with their ranking in the daily 
subnetworks. We did so as follows: for a given underlay, 
we ranked each node according to its degree, with the 
most highly connected node being given rank 0; nodes 
with the same degree were grouped together and given 
a rank which was the average of the range they covered. 
We then generated 10 daily subnetworks of S=l million 
steps each, and averaged the rank each node had over 
all ten subnetworks, dealing with ties in the same way 
as above. We did this for C = and C — 10, and 
Fig. [5] shows the underlay rank plotted versus the average 
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daily rank so calculated. When multiple nodes have the 
same underlay rank, the average daily rank has itself been 
averaged over all nodes with the same underlay rank. In 
the C — case, we see that the underlay and daily ranks 
have a strong positive correlation: the daily rank of a 
node always depends strongly on its rank in the underlay. 
The C = 10 curve is much different: the average daily 
rank is roughly constant, once we consider nodes below 
the top 50,000 (5%) or so. In addition to the average 
daily rank, we can also consider the best (i.e. minimum) 
daily rank that a node reaches over the 10 daily networks; 
in doing so, we see (in the inset to Fig. [5]) that every node 
in the C > case has an opportunity to be among the 
top 1500 nodes (0.15%) in some daily network; without 
preferential attachment, most nodes can't crack the top 
10,000 (1%). The C > case is extremely egalitarian 
in this regard: any node, even those with the smallest 
degree in the underlay, can be among the most important 
nodes in some daily network. This behavior agrees with 
the empirical results reported in [33] 23J . 



this dissimilarity as a function of the aggregate size M. 
For the pure random walk case (C = 0), the dissimilarity 
decreases exponentially (~ e -M/M ^ w ^ m q ~ —7.5), 
so that after 20 days the top nodes in both networks are 
basically identical. In contrast, the aggregate generated 
with a strong dynamic preferential attachment (C = 10) 
has a dissimilarity which decays much more slowly, and 
possibly without a characteristic timescale. This is the 
behavior seen in empirical DSF networks [12 03], and 
hints at a "multiscale" structure for the network dy- 
namics, where networks at each timescale form scale-free 
topologies while the specific links in existence vary dra- 
matically between observation timescales as well as over 
time. To alleviate concerns that this disparity may be a 
finite-size effect, we include data in Fig. [6] for networks 
of N = 50,000 nodes and N = 2,000,000 nodes (where 
S = N in all cases): increasing the network size does 
nothing to change the C = case, and actually results 
in an even slower rate of decay for the reinforced C = 10 
networks. 
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FIG. 6. Dissimilarity of networks aggregated over times rang- 
ing from 1 to 50 days. Dissimilarity is defined as the Top-2% 
overlap O (Eq. 2) subtracted from 100%, and is shown for 
both pure random walks (C = 0) and strongly weighted walks 
(C = 10). We show results for networks with 50 thousand, 1 
million, and 2 million nodes to show how the system scales; 
all walks have as many steps as there are nodes (S = TV). 

One last empirical result seen in [42] involves the ag- 
gregate of M daily subnetworks, which are formed by 
connecting two nodes together if they are connected in 
any of the M subnetworks. Given their independence 
and random starting points, the aggregate of M daily 
networks must eventually converge to the underlay as M 
grows large, but whether this convergence occurs over 
a particular timescale Mo, or whether it converges in a 
scale-free manner, is an interesting question. One way 
to measure this convergence is by calculating the Top- 
2% overlap O between the aggregate and the underlay, 
and subtracting this from 100%; if this quantity (called 
the "dissimilarity" reaches zero, then both the aggregate 
and the underlay both have the same set of hubs, and are 
roughly (though not precisely) equivalent. Fig. [6] shows 



IV. CONCLUSION 

We have embarked on a research program designed 
to develop universal models that can recreate empiri- 
cally observed phenomena in dynamic complex networks. 
We have shown that, using a suitable reinforced random 
walk on a "long-term" underlay network, one is able to 
produce instantaneous networks which reproduce qual- 
itatively characteristic features of real world dynamic 
networks. This includes, in particular, the construc- 
tion of scale-free sub-networks of a scale-free "underlay" 
network, whose local hubs substantially differ from sub- 
network to sub-network and from those of the underlay. 
We have presented evidence that dynamic preferential 
attachment (as opposed to a pure random walk) is a nec- 
essary requirement of dynamic scale-free networks: the 
preferential attachment mechanism limits the walk's ex- 
ploration of the network, which gives early-visited nodes 
higher ranks than they would have with a more random 
traversal. We hope our model will stimulate further em- 
pirical and theoretical work, and provide a framework 
for analyzing the influence of link/structural dynamics 
on dynamical processes on complex networks. 

Appendix A: Calculating the top— n overlap 

If there are no ties for nth place in the network, then 
the top-n overlap of subnetworks A and B is given by 
the formula 

o = - 1 n B top \ , (Al) 
n 

where A top is a set of the top n nodes in network A, as 
Stop is for 23, and |-| denotes the number of elements in a 
set. Given the finite size of our networks, however, there 
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are generally large groups of nodes with the same degree, 
making the choice of the top n nodes somewhat arbitrary. 
To account for this, we modify Eq. |A1| to prorate those 
nodes which are tied for the nth slot, as follows: let a 
and b be the degree of the nth highest node in networks 
A and B, respectively. We define the set Ay to be the 
set of nodes in A with degree greater than a, while A— 
is the (guaranteed nonempty) set of nodes in A with de- 
gree equal to a; we define By and B = similarly. Imagine 
constructing the list A top of the top n nodes: it will auto- 
matically contain the entire set Ay, and will have room 
for n — \A> \ of the nodes in A = . If we choose nodes at 
random from A = to fill those available slots, then each 



node from A = has probability 

n-|A>| 



A 



\AJ 



(A2) 



of appearing in the list A top . (The same applies for net- 
work B.) Thus, when calculating the overlap between 
A and B, we will count each node in A = and B = as 
being worth a fraction of a regular node — fj± and /s, 
respectively — define the modified overlap 0(A, B) to be 



O 



l r 



|A > ns>| + / A |^ = nB > | 



+f B \AyC\B = \ + f A f B \A = ^B = \ (A3) 



Note that if there is no tie for nth place in either 
network, then !a = !b = 1 , A t op = (Ay U A = ), and 



B 



top 



Eq. Al described above. 



(By U B = ); and Eq. A3 reduces to the simpler 
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